QTL Mapping Using RIL Population.

Genetic maps are an excellent tool for the analysis of important traits, the development of which is the result of the combined expression of several genes, enabling the genomic localization of the factors determining them. Such features, characterized by a normal distribution of values, are referred to as quantitative or polygenic. The analysis of their genetic background using a chromosome map is called the mapping of quantitative traits loci (QTL). QTL analysis is a statistical method of determining the genetic association of phenotypic data (trait measurements) with genotypic data (DNA markers assigned to linkage groups).There are numerous tools developed for QTL mapping. This chapter introduces Windows QTL Cartographer with Composite Interval Mapping (CIM) method, which estimates the QTL position by combining interval mapping with multiple regression. The genotypic and phenotypic data used in the exemplary QTL mapping procedure were obtained for the recombinant inbred line (RIL) population of rye. Plant height, assessed in three seasons, was the exemplary trait under study.

Genetic Map Construction Using F2 and RIL/DH Mapping Population.

Genetic mapping is the determination of the position and relative genetic distance between genes or molecular markers in the chromosomes of a particular species. The construction of genetic maps uses data from the genotyping of the mapping population. Among the different mapping populations used, two are relatively common: the F2 and recombinant inbred lines (RILs) obtained as a result of the controlled crossing of genetically diverse parental forms (e.g., inbred lines). Also, the dihaploid (DH) population is often used in plants, but obtaining DHs in different crops, including rye, is very difficult or even impossible. Any molecular marker system can be used for genotyping. Polymorphic markers are used for linkage analysis, differentiating parental forms with segregation in the mapping population, consistent with the appropriate single-gene model. A genetic map is a great source of information on a species and can be an exquisite tool for analyzing important quantitative traits (QT).This chapter presents the procedure of genetic map construction with two different algorithms using the JoinMap5.0 program. First, the Materials section briefly informs about the mapping program, showing how to obtain a mapping population and prepare data for mapping. Finally, the Methods section describes the protocol for the mapping procedure itself.

In search of the relationship between the rye polyamine oxidase (PAO) gene and resistance to powdery mildew (PM).

Powdery mildew (PM), a common cereal disease in cultivated areas, including Europe and other temperate regions, is caused by the fungus Blumeria graminis. While PM is one of the most important wheat leaf diseases globally, rye is highly tolerant to PM. It has been reported that in barley infected with PM, polyamine oxidase (PAO) activity related to the production of hydrogen peroxide (H2O2) has increased, which may promote defense against biotrophic or hemibiotrophic pathogens. The current study aimed to assess the relationship between the segregation of the polymorphic marker for rye PAO (ScPAO) and the level of PM infection in plants. The genetic mapping in two interline populations shows that ScPAO is located on chromosome 7R. Further analysis comparing ScPAO location to mapped wheat (Triticum aestivum L.) PAO duplicates suggests the ScPAO homology with TaPAO6 or TaPAO7. A possible association of ScPAO from 7R with PM resistance is demonstrated in the recombinant inbred lines (RIL)-L population phenotyped for PM infection. Finally, three novel QTLs for PM resistance on the 7R chromosome of rye are detected.

Variation between glaucous and non-glaucous near-isogenic lines of rye (Secale cereale L.) under drought stress.

Glaucous (811, L35, and RXL10) and non-glaucous (811bw, L35bw, and RXL10bw) near-isogenic lines (NILs) of rye (Secale cereale L.) forming three pairs of inbred lines were the subject of the research. The research aimed to study the relationship between wax cover attributes and the physio-biochemical drought reactions and yield of rye NILs and to uncover the differences in drought resistance levels of these lines. The greatest differences between glaucous and non-glaucous NILs were observed in the RXL10/RXL10bw pair. Of particular note were the stable grain number and the thousand grain weight of the non-glaucous line RXL10bw under drought and the accompanying reactions, such as an approximately 60% increase in MDA and a two-fold increase in wax amount, both of which were significantly higher than in the glaucous line RXL10 and in other NILs. The surprisingly high level of MDA in the RXL10bw line requires further analysis. Moreover, additional wax crystal aggregates were found under drought conditions on the abaxial leaf surface of the glaucous lines 811 and RXL10. The use of rye NILs indicated that line-specific drought resistance could be associated with wax biosynthetic pathways involved in physiological and biochemical responses important for increased drought resistance.

Functional Genomics for Plant Breeding.

To face the rapidly growing world human population, an increase in agricultural productivity and production is necessary to overcome the enhanced food demand [...].

Correlation between leaf epicuticular wax composition and structure, physio-biochemical traits and drought resistance in glaucous and non-glaucous near-isogenic lines of rye.

The objective of this research was to investigate the differences between glaucous and non-glaucous near-isogenic lines (NILs) of winter rye (Secale cereale L.) in terms of epicuticular wax layer properties (weight, composition, and crystal morphology), selected physiological and biochemical responses, yield components, above-ground biomass, and plant height under soil drought stress. An important aspect of this analysis was to examine the correlation between the above characteristics. Two different NIL pairs were tested, each consisting of a typical glaucous line and a non-glaucous line with a recessive mutation. The drought experiment was conducted twice (2015-2016). Our study showed that wax accumulation during drought was not correlated with higher leaf hydration and glaucousness. Environmental factors had a large impact on the response of the lines to drought in individual years, both in terms of physiological and biochemical reactions, and the composition of epicuticular leaf wax. The analysed pairs displayed significantly different responses to drought. Demonstration of the correlation between the components of rye leaf wax and the physiological and biochemical parameters of rye NILs is a significant achievement of this work. Interestingly, the study showed a correlation between the wax components and the content of photosynthetic pigments and tocopherols, whose biosynthesis, similarly to the biosynthesis of wax precursors, is mainly located in chloroplasts. This suggests a relationship between wax biosynthesis and plant response to various environmental conditions and drought stress.

Chromosome-scale genome assembly provides insights into rye biology, evolution and agronomic potential.

Rye (Secale cereale L.) is an exceptionally climate-resilient cereal crop, used extensively to produce improved wheat varieties via introgressive hybridization and possessing the entire repertoire of genes necessary to enable hybrid breeding. Rye is allogamous and only recently domesticated, thus giving cultivated ryes access to a diverse and exploitable wild gene pool. To further enhance the agronomic potential of rye, we produced a chromosome-scale annotated assembly of the 7.9-gigabase rye genome and extensively validated its quality by using a suite of molecular genetic resources. We demonstrate applications of this resource with a broad range of investigations. We present findings on cultivated rye's incomplete genetic isolation from wild relatives, mechanisms of genome structural evolution, pathogen resistance, low-temperature tolerance, fertility control systems for hybrid breeding and the yield benefits of rye-wheat introgressions.

How Machine Learning Methods Helped Find Putative Rye Wax Genes Among GBS Data.

The standard approach to genetic mapping was supplemented by machine learning (ML) to establish the location of the rye gene associated with epicuticular wax formation (glaucous phenotype). Over 180 plants of the biparental F2 population were genotyped with the DArTseq (sequencing-based diversity array technology). A maximum likelihood (MLH) algorithm (JoinMap 5.0) and three ML algorithms: logistic regression (LR), random forest and extreme gradient boosted trees (XGBoost), were used to select markers closely linked to the gene encoding wax layer. The allele conditioning the nonglaucous appearance of plants, derived from the cultivar Karlikovaja Zelenostebelnaja, was mapped at the chromosome 2R, which is the first report on this localization. The DNA sequence of DArT-Silico 3585843, closely linked to wax segregation detected by using ML methods, was indicated as one of the candidates controlling the studied trait. The putative gene encodes the ABCG11 transporter.

The GAMYB gene in rye: sequence, polymorphisms, map location, allele-specific markers, and relationship with α-amylase activity.

BACKGROUND: Transcription factor (TF) GAMYB, belonging to MYB family (named after the gene of the avian myeloblastosis virus) is a master gibberellin (GA)-induced regulatory protein that is crucial for development and germination of cereal grain and involved in anther formation. It activates many genes including high-molecular-weight glutenin and α-amylase gene families. This study presents the first attempt to characterize the rye gene encoding GAMYB in relation to its sequence, polymorphisms, and phenotypic effects. RESULTS: ScGAMYB was mapped on rye chromosome 3R using high-density Diversity Arrays Technology (DArT)/DArTseq-based maps developed in three mapping populations. The ScGAMYB sequences were identified in RNA-seq libraries of four rye inbred lines. The transcriptome used for the search contained almost 151,000 transcripts with a median contig length of 500 nt. The average amount of total base raw data was approximately 9 GB. Comparative analysis of the ScGAMYB sequence revealed its high level of homology to wheat and barley orthologues. Single nucleotide polymorphisms (SNPs) detected among rye inbred lines allowed the development of allele specific-PCR (AS-PCR) markers for ScGAMYB that might be used to detect this gene in wide genetic stocks of rye and triticale. Segregation of the ScGAMYB alleles showed significant relationship with α-amylase activity (AMY). CONCLUSIONS: The research showed the strong similarity of rye GAMYB sequence to its orthologues in other Graminae and confirmed the position in the genome consistent with the collinearity rule of cereal genomes. Concurrently, the ScGAMYB coding sequence (cds) showed stronger variability (24 SNPs) compared to the analogous region of wheat (5 SNPs) and barley (7 SNPs). The moderate regulatory effect of ScGAMYB on AMY was confirmed, therefore, ScGAMYB was identified as a candidate gene for partial control of α-amylase production in rye grain. The predicted structural protein change in the aa region 362-372, caused by a single SNP (C/G) at the 1100 position in ScGAMYB cds and single aa sequence change (S/C) at the 367 position, is the likely cause of the differences in the effectiveness of ScGAMYB regulatory function associated with AMY. The development of sequence-based, allele-specific (AS) PCR markers could be useful in research and application.

Genetic Parameters and QTLs for Total Phenolic Content and Yield of Wheat Mapping Population of CSDH Lines under Drought Stress.

A doubled haploid population of 94 lines from the Chinese Spring × SQ1 wheat cross (CSDH) was used to evaluate additive and epistatic gene action effects on total phenolic content, grain yield of the main stem, grain number per plant, thousand grain weight, and dry weight per plant at harvest based on phenotypic and genotypic observations of CSDH lines. These traits were evaluated under moderate and severe drought stress and compared with well-watered plants. Plants were grown in pots in an open-sided greenhouse. Genetic parameters, such as additive and epistatic effects, affecting total phenolic content, were estimated for eight year-by-drought combinations. Twenty-one markers showed a significant additive effect on total phenolic content in all eight year-by-drought combinations. These markers were located on chromosomes: 1A, 1B, 2A, 2B, 2D, 3A, 3B, 3D, 4A, and 4D. A region on 4AL with a stable QTL controlling the phenolic content, confirmed by various statistical methods is particularly noteworthy. In all years and treatments, three markers significantly linked to QTLs have been identified for both phenols and yield. Thirteen markers were coincident with candidate genes. Our results indicated the importance of both additive and epistatic gene effects on total phenolic content in eight year-by-drought combinations.

Putative candidate genes responsible for leaf rolling in rye (Secale cereale L.).

BACKGROUND: Rolling of leaves (RL) is a phenomenon commonly found in grasses. Morphology of the leaf is an important agronomic trait in field crops especially in rice; therefore, majority of the rice breeders are interested in RL. There are only few studies with respect to RL of wheat and barley; however, the information regarding the genetic base of RL with respect to the shape of leaf in rye is lacking. To the best of our knowledge, this is the first study on the localization of loci controlling RL on high density consensus genetic map of rye. RESULTS: Genotypic analysis led to the identification of 43 quantitative trait loci (QTLs) for RL, grouped into 28 intervals, which confirms the multigenic base of the trait stated for wheat and rice. Four stable QTLs were located on chromosomes 3R, 5R, and 7R. Co-localization of QTL for RL and for different morphological, biochemical and physiological traits may suggests pleiotropic effects of some QTLs. QTLs for RL were associated with QTLs for such morphological traits as: grain number and weight, spike number per plant, compactness of spike, and plant height. Two QTLs for RL were found to coincide with QTLs for drought tolerance (4R, 7R), two with QTLs for heading earliness (2R, 7R), one with α-amylase activity QTL (7R) and three for pre-harvest sprouting QTL (1R, 4R, 7R). The set of molecular markers strongly linked to RL was selected, and the putative candidate genes controlling the process of RL were identified. Twelve QTLs are considered as linked to candidate genes on the base of DArT sequences alignment, which is a new information for rye. CONCLUSIONS: Our results expand the knowledge about the network of QTLs for different morphological, biochemical and physiological traits and can be a starting point to studies on particular genes controlling RL and other important agronomic traits (yield, earliness, pre-harvest sprouting, reaction to water deficit) and to appoint markers useful in marker assisted selection (MAS). A better knowledge of the rye genome and genes could both facilitate rye improvement itself and increase the efficiency of utilizing rye genes in wheat breeding.

Effective BAC clone anchoring with genotyping-by-sequencing and Diversity Arrays Technology in a large genome cereal rye.

Identification of bacterial artificial chromosome (BAC) clones containing specific sequences is a prerequisite for many applications, such as physical map anchoring or gene cloning. Existing BAC library screening strategies are either low-throughput or require a considerable initial input of resources for platform establishment. We describe a high-throughput, reliable, and cost-effective BAC library screening approach deploying genotyping platforms which are independent from the availability of sequence information: a genotyping-by-sequencing (GBS) method DArTSeq and the microarray-based Diversity Arrays Technology (DArT). The performance of these methods was tested in a very large and complex rye genome. The DArTseq approach delivered superior results: a several fold higher efficiency of addressing genetic markers to BAC clones and anchoring of BAC clones to genetic map and also a higher reliability. Considering the sequence independence of the platform, the DArTseq-based library screening can be proposed as an attractive method to speed up genomics research in resource poor species.

Bidirectional selective genotyping approach for the identification of quantitative trait loci controlling earliness per se in winter rye (Secale cereale L.).

The genes controlling earliness of plants include genes responsible for vernalisation (Vrn) and photoperiod (Ppd), and those that are not entirely associated with a response to temperature or light. The last group of loci is known as earliness per se (Eps). Eps genes have been most commonly reported in the scientific literature as quantitative trait loci (QTL). The objective of this study was to use a bidirectional selective genotyping (BSG) method for the identification of loci controlling Eps, assessed at the heading stage in a newly developed rye population of recombinant inbred lines. We identified four linkage groups representing fragments of chromosomes, 1R, 5R, 6R and 7R, and containing, in total, 155 DArT markers significantly associated with earliness in rye. The results are discussed in the context of previously published QTL analyses.

Genetic mapping of the ScHd1 gene in rye and an assessment of its relationship with earliness per se and plant morphology.

A fragment of the ScHd1 gene derived from eight inbred lines was sequenced and showed homology to other Hd1 genes from different cereals. Sequences were analysed with respect to the presence of a single-nucleotide polymorphism (SNP) difference. A C-T transition at position 312 of the consensus sequence was found, which distinguished two lines from the remaining six. The deduced amino acid sequence revealed a high identity (93%) to a Hd1-like protein from wheat. The identified mutation allowed the localisation of ScHd1 on a genetic map of rye (6RS). A small, statistically significant linkage between ScHd1 and earliness per se (eps) and some morphological traits was also established. The chromosomal region, including the S76 allele for the ScHd1 gene was linked to earlier heading, elongated spikes, a greater number of spikelets per spike and an increased weight of 1000 kernels.

The application of high-density genetic maps of rye for the detection of QTLs controlling morphological traits.

The development of genetic maps is, nowadays, one of the most intensive research activities of plant geneticists. One of the major goals of genome mapping is the localisation of quantitative trait loci (QTLs). This study was aimed at the identification of QTLs controlling morphological traits of rye and comparison of their localisation on genetic maps constructed with the use of genetically different germplasms. For QTL analyses, two high-density consensus maps of two populations (RIL-S and RIL-M) of recombinant inbred lines (RIL) were applied. Plant height (Ph), length of spikes (Sl) and the number of spikelets per spike (Sps) were studied in both populations. Additionally, the number of kernels per spike under isolation (Kps), the weight of kernels per spike (Kw) and thousand kernel weight (Tkw) were assessed in the RIL-M population. Except for Tkw, the majority of the traits were correlated to each other. The non-parametric Kruskal-Wallis (K-W) test and composite interval mapping (CIM) revealed 18/48 and 24/18 regions of rye chromosomes engaged in the determination of Ph, Sl and Sps in the RIL-S and RIL-M populations, respectively. An additional 18/15 QTLs controlling Kps, Kw and Tkw were detected on a map of the RIL-M population. A numerous group of QTLs detected via CIM remained in agreement with the genomic regions found when the K-W test was applied. Frequently, the intervals indicated by CIM were narrower.

The importance of chromosomes from the sixth homeologic group in the restoration of male fertility in winter triticale with Triticum timopheevii cytoplasm.

The sterilising cytoplasm from Triticum timopheevii is presently considered to be the most promising as regards to the seed production of triticale hybrid cultivars. This study was aimed at the utilisation of Diversity Arrays Technology (DArT) for the preliminary identification of genomic regions with loci controlling male sterility/fertility in the field-grown F2 generation of the interline hybrid between male sterile line CMS-Salvo 15/1 and restorer line Stan I. The fertility of plants was examined by visual scoring as well as by the assessment of seed setting within bagged spikes. For DNA analyses, 92 individuals representing opposite phenotypes (male sterile vs. fully male fertile) were chosen from the whole F2 population, which consisted of 414 plants. The constructed genetic map consists of 759 DArT markers distributed in 24 linkage groups that cover a distance of 974.4 cM. Application of the interval mapping method and the Kruskal-Wallis test enabled the identification of six genomic regions engaged in the restoration of male fertility within the mapping population. The most effective restorer genes were found on chromosomes of the sixth homeologic group, i.e. on 6R (the most efficient), 6A and 6B. Additionally, linkage groups assigned to chromosomes 1BS, 3A and 3A/3B were important for the determination of male fertility.

Detection of the quantitative trait loci for α-amylase activity on a high-density genetic map of rye and comparison of their localization to loci controlling preharvest sprouting and earliness.

The objectives of the research were to determine the position of quantitative trait loci (QTL) for α-amylase activity on the genetic map of a rye recombinant inbred line population-S120 × S76-and to compare them to known QTL for preharvest sprouting and heading earliness. Fourteen QTL for α-amylase activity on all seven chromosomes were identified. The detected QTL were responsible for 6.09-23.32% of α-amylase activity variation. The lowest LOD value (2.22) was achieved by locus QAa4R-M3 and the highest (7.79) by locus QAa7R-M1. Some QTL intervals for features of interest overlapped partially or completely. There were six overlapping QTL for α-amylase activity and preharvest sprouting (on 1R, 3R, 4R, 6R, 7R) and the same number for preharvest sprouting and heading earliness (on 1R, 2R, 6R, 7R). Furthermore, there was one interval partially common to all three traits, mapped on the long arm of chromosome 1R. Testing of lines originating from hybrid breeding programs, such as S120 and S76, may provide important information about the most significant genes and markers for selection in commercial breeding. Among the statistically significant markers selected in the Kruskal-Wallis test (P < 0.005), there were 55 common ones for preharvest sprouting and heading earliness (1R, 2R, 6R), 30 markers coinciding between α-amylase activity and preharvest sprouting (5R, 7R) and one marker for α-amylase activity and heading earliness (6R). ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-011-9627-1) contains supplementary material, which is available to authorized users.

A high density consensus map of rye (Secale cereale L.) based on DArT markers.

BACKGROUND: Rye (Secale cereale L.) is an economically important crop, exhibiting unique features such as outstanding resistance to biotic and abiotic stresses and high nutrient use efficiency. This species presents a challenge to geneticists and breeders due to its large genome containing a high proportion of repetitive sequences, self incompatibility, severe inbreeding depression and tissue culture recalcitrance. The genomic resources currently available for rye are underdeveloped in comparison with other crops of similar economic importance. The aim of this study was to create a highly saturated, multilocus linkage map of rye via consensus mapping, based on Diversity Arrays Technology (DArT) markers. METHODOLOGY/PRINCIPAL FINDINGS: Recombinant inbred lines (RILs) from 5 populations (564 in total) were genotyped using DArT markers and subjected to linkage analysis using Join Map 4.0 and Multipoint Consensus 2.2 software. A consensus map was constructed using a total of 9703 segregating markers. The average chromosome map length ranged from 199.9 cM (2R) to 251.4 cM (4R) and the average map density was 1.1 cM. The integrated map comprised 4048 loci with the number of markers per chromosome ranging from 454 for 7R to 805 for 4R. In comparison with previously published studies on rye, this represents an eight-fold increase in the number of loci placed on a consensus map and a more than two-fold increase in the number of genetically mapped DArT markers. CONCLUSIONS/SIGNIFICANCE: Through the careful choice of marker type, mapping populations and the use of software packages implementing powerful algorithms for map order optimization, we produced a valuable resource for rye and triticale genomics and breeding, which provides an excellent starting point for more in-depth studies on rye genome organization.

DArT markers tightly linked with the Rfc1 gene controlling restoration of male fertility in the CMS-C system in cultivated rye (Secale cereale L.).

The Rfc1 gene controls restoration of male fertility in rye (Secale cereale L.) with sterility-inducing cytoplasm CMS-C. Two populations of recombinant inbred lines (RIL) were used in this study to identify DArT markers located on the 4RL chromosome, in the close vicinity of the Rfc1 gene. In the population developed from the 541×2020LM intercross, numerous markers tightly linked with the restorer gene were identified. This group contained 91 DArT markers and three SCARs additionally analyzed in the study. All these markers were mapped in the distance not exceeding 6 cM from the gene of interest. In the second mapping population (541×Ot1-3 intercross), only 9 DArT markers located closely to the Rfc1 gene were identified. Five of these DArT markers were polymorphic in both populations.

A consensus map of chromosome 6R in rye (Secale cereale L.).

Four F(2) mapping populations derived from crosses between rye inbred lines DS2 x RXL10, 541 x Ot1-3, S120 x S76 and 544 x Ot0-20 were used to develop a consensus map of chromosome 6R. Thirteen marker loci that were polymorphic in more than one mapping population constituted the basis for the alignment of the four maps using the JoinMap v. 3.0 software package. The consensus map consists of 104 molecular marker loci including RFLPs, RAPDs, AFLPs, SSRs, ISSRs, SCARs, STSs and isozymes. The average distance between the marker loci is 1.3 cM, and the total map length is 135.5 cM. This consensus map may be used as a source of molecular markers for the rapid development of new maps of chromosome 6R in any mapping population.